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This study is part of a research program aimed to investigate the relations between instantons, 
monopoles, and chiral symmetry breaking. Monopoles are important 3-dimensional topological 
configurations existing in QCD, which are believed to produce colour confinement. Instantons are 
4-dimensional topological configurations and are known to be related to chiral symmetry breaking. 
To study the relation between monopoles and instantons we generate configurations adding to the 
vacuum state static monopole-antimonopole pairs of opposite charges by use of a monopole creation 
operator. We observe that the monopole creation operator only adds long monopole loops to the 
configurations. We then count the number of fermion zero modes using Overlap fermions as a tool. 
As a result we find that each monopole-antimonopole pair of magnetic charge one adds one zero 
mode of chirality ±1, i.e. one instanton of topological charge ±1. 

PACS numbers: 11.15.Ha, 12.38.Gc 


I. INTRODUCTION 

Monopoles are 3-dimensional topological QCD config¬ 
urations, which are most probably responsible for colour 
confinement, by condensing in the vacuum and producing 
dual superconductivity [l|. This mechanism has exten¬ 
sively been studied on the Lattice during the years with 
different approaches. 

Instantons are 4-dimensional configurations and are re¬ 
lated to Chiral symmetry breaking as described e.g. by 
the instanton liquid model Q. 

There exist strong hints that the role of the two con¬ 
figurations can be related. In Ref. Q it is shown that 
monopoles catalyse instantons, in the sense that the 
probability of having an instanton configuration in the 
presence of a monopole is of order 1. The analysis 
there was specifically addressed to grand-unified mod¬ 
els in view of the proton decay, but it is equally valid in 
QCD. On the other hand the work of Ref. ai shows 
that instantons, and more generally calorons are made of 
monopoles. 

As a first step towards a quantitative understanding of 
this relation, we add static monopole-antimoncmole pairs 
to the vacuum of an SU{S) pure gauge theory @-@1, with 
the aim of changing the number of monopoles without 
significant modifications of the vacuum state. Indeed, 
if the vacuum is a dual superconductor adding a pair 
of magnetic charges does not modify the ground state, 
which is an eigenstate of the creation operators, except 
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in the total number of positive and negative monopoles. 
Magnetic charges are in any case shielded and the addi¬ 
tional monopoles are one-dimensional structures prop¬ 
agating in time: their effect on the vacuum vanishes 
at large volumes as V~^. Moreover the total mag¬ 
netic charge being zero the operator which creates the 
monopoles has no infrared problems Q. On the states 
constructed in this way we measure the change in the 
number of instantons. 

As a tool to count instantons we use the zero modes of 
an Overlap Dirac Operator. In doing that we realise that 
each lattice configuration contains zero modes of only one 
sign, either -|- or —, and never zero modes of opposite sign 
at the same time: this contrasts with the naive expecta¬ 
tion that, zero modes being localised and the correlation 
length finite, independent parts of the lattice could con¬ 
tain instantons of opposite signs. We also would expect 
that the number of instantons and of anti-instantons, i.e. 
the number of zero modes, be proportional to the volume. 
We argue that, for some reason to be quantitatively stud¬ 
ied and understood we only observe the net topological 
charge, while pairs of zero modes of opposite sign escape 
detection. Similar views can be found in Q. Maybe this 
phenomenon is due to our rather small volumes, and will 
disappear at larger volumes. We check this argument by 
doing a careful study of the distributions of topological 
charges, a careful determination of the topological sus¬ 
ceptibility, which we compare with existing results in the 
literature, to convince ourselves that the determination 
of the topological charge is correct and also that our code 
is correct [Section II]. 

The number of instantons can in any case be deter¬ 
mined from the topological susceptibility under the very 
general assumptions of translation invariance, existence 
of a finite correlation length, and CP invariance [Section 
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III], 

In Section IV we briefly recall the definition of the 
operator which creates a static pair of monopoles propa¬ 
gating in time from t = —oo to t = -|-oo and sitting at a 
given spatial distance. 

In Section V we check that the operator really cre¬ 
ates monopoles by using the technique of Ref. [ifl [ll| . 
We find that indeed it creates so called long monopoles, 
which are not short-range fluctuations but wrap the lat¬ 
tice along the time direction. 

Finally in Section VI we look at the configurations with 
monopole pairs of different magnetic charges added, we 
measure the topological charge and from it we compute 
the number of instantons. The result is that a pair of 
monopoles of charges ±1 produces one instanton or anti- 
instanton. A pair with magnetic charge rric produces rric 
instantons or anti-instantons. 

In Section VII we draw some conclusions. 

Some of the results have already been reported at the 
LATTICE 2014 Conference [l^ and at the Conference 
Confinement 2014 [l3l| . 


A. Lattice spacing 


To fix the scale, we determine the lattice spacing by 
using of the analytic interpolation of Ref. m. As a check 
we also measure the string tension and the qq potential 
from Wilson loops. We use APE smearing 18| of the 
spacial link variables to suppress excited states and fit 
the function V{R) = a ■ R — a/R + C to the potential 
V(R ). The distance R between quark and anti-quark is 
improved to Ri using the Green function as defined in 
Ref. [13111. The Jackknife method is used to estimate 
the statistical errors. 

The lattice spacing is then computed in two different 
ways (1) from a and Sommer scale tq = 0.5 [fm], 
and (2) from the String tension y/a = 440 [MeV] 
[Table H]. The results are reasonably consistent with the 
analytic interpolation. Therefore, we take the lattice 
spacing from the analytic interpolation [17|, and use the 
Sommer scale Tq = 0.5 [fm]. 


B. Simulation details 


II. OVERLAP FERMIONS 

Wilson fermions explicitly break Chiral symmetry. 
Therefore we use Overlap fermions which preserve Chi¬ 
ral symmetry on the Lattice. The Ginsparg-Wilson re¬ 
lation [l3[ describes Chiral symmetry in Lattice gauge 
theory, 

75 II -I- £>75 = aDRj^D. (1) 

a is the lattice spacing, D the Overlap Dirac operator, 
and R a parameter. The right hand side of this relation 
is non zero, because of the Nielsen-Ninomiya [ll no go 
theorem. Multiplying both sides of Eq. (HD by the inverse 
of the Overlap Dirac operator, D~^, gives 


We generate configurations using the Wilson gauge ac¬ 
tion and periodic boundary conditions. The number of 
iterations for the thermalization is 0 ( 2.0 x 10 "'^), and the 
configurations are sampled after 0(5.0 x 10^) iterations 
between them. The simulation parameters are listed in 
Table mi 

We construct the Overlap Dirac operator D from the 
massless Wilson fermions Dw which are computed with 
the gauge links of the configurations. An anti-periodic 
boundary condition in the temporal direction is used for 
the Wilson fermions. 

In more detail, we carry out the numerical compu¬ 
tations by the technique of Ref. [2fll - [^ . The massless 
Overlap Dirac operator D{p) is 


75 D ^+D i75=ai?75, ( 2 ) 

showing that Chiral symmetry breaking of the propaga¬ 
tor D~^ is an operator of 0{a) vanishing in the contin¬ 
uum limit. The exact form of the Overlap Dirac opera¬ 
tor is defined in Ref. (T^j in terms of the massless Wilson 
Dirac operator Dw- 

Ra V y/WAj 

A = —Mq -t- aDw, Mq is a parameter, 0 < Mq < 2. Dw 
is the massless Wilson Dirac operator defined as follows: 

Dw = l{MK + V^)-aV;y,} (4) 

[V*V’](?^) = ^ - /i)| (6) 


D{p) 


P 

a 



Dwjp) 

\J Dw{p)^Dw{p) 


( 7 ) 


Dw ip) is computed from massless Wilson Dirac operator 
Dw, Eq- O, as follows: 

Dw{p) = Dw — ^ ( 8 ) 

p is a (negative) mass parameter, 0 < p < 2. We set 
p = 1.4 in this study. Using the Hermitian Wilson Dirac 
operator Hwip) = j^Dwip), the sign function is defined 
as in Ref. M 


e{Hw{p)) = 


Hwjp) 

\/Hw{pVHw{p) 


( 9 ) 


Therefore, the massless Overlap Dirac operator is com¬ 
puted by the sign function 

D{p) = ^{l + 'l5£{Hw{p))} ■ (10) 
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TABLE I. Determination of the lattice spacing. The lattice spacing is computed by the analytic function of [O]- The 
lattice spacing is computed by our simulations of the qq potential (n is the number of smearing steps and a is the weight 
factor of smearing, FR is the range of the fit) and from our determination of the string tension. 


p 

[fm] 

[fm] 

a*-2) [fjji] 

V 

(n, a) 

FR {Ri/a) 

X^/ndf 

N f 

^ ’ conj. 

6.00 

9.315 X 10-2 

9.39(9) X 10-2 

9.93(8) X 10-2 

18'‘ 

(20, 0.5) 

0.9 - 7.0 

4.3/4.0 

440 


TABLE 11. This table summarises the results 

/ndf refers to the ht of the function Pq = 

of our simulation; 

2i 

s for the number of zero modes and the topological susceptibility. 

listributions of the topological charges. 


P 

a/ro 

V 

VIA 

Zero 


{QAAIv 

X^/ndf 

TV f 

5.79 

0.2795 

12"^ 

126.5 

2.32(12) 

8.3(7) 

0.066(6) 

4.7/14.0 

200 

5.81 

0.2659 

10"^ 

49.96 

1.44(4) 

3.50(17) 

0.070(3) 

8.5/11.0 

844 



14" 

191.9 

2.80(13) 

12.2(1.1) 

0.064(5) 

12.4/17.0 

249 



16" 

327.4 

3.63(17) 

20.9(1.8) 

0.064(6) 

19.6/19.0 

275 

5.85 

0.2484 

12" 

78.95 

1.87(10) 

5.3(5) 

0.068(6) 

3.9/11.0 

200 



16" 

249.5 

3.30(13) 

17.0(1.3) 

0.068(5) 

22.0/19.0 

344 

5.86 

0.2395 

14" 

126.5 

2.32(9) 

8.3(7) 

0.066(6) 

16.7/15.0 

338 

5.90 

0.2216 

12" 

49.96 

1.40(4) 

3.32(16) 

0.067(3) 

10.1/11.0 

835 



16" 

157.9 

2.48(11) 

10.3(9) 

0.065(5) 

18.1/19.0 

320 

5.93 

0.2129 

14" 

78.95 

1.95(9) 

5.8(5) 

0.073(6) 

15.0/12.0 

278 

5.99 

0.1899 

14" 

49.96 

1.39(4) 

3.17(15) 

0.063(3) 

8.8/11.0 

862 

6.00 

0.1863 

12" 

24.98 

0.88(4) 

1.47(11) 

0.059(4) 

6.9/8.0 

430 



14" 

46.28 

1.35(5) 

3.17(19) 

0.069(4) 

10.2/11.0 

592 



12^ X 24 

49.96 

1.34(4) 

2.98(15) 

0.060(3) 

11.7/12.0 

790 



16" 

78.95 

1.81(7) 

5.4(4) 

0.068(5) 

7.0/14.0 

405 



18" 

126.5 

2.42(11) 

8.8(7) 

0.069(5) 

16.7/14.0 

258 

6.07 

0.1662 

16" 

49.96 

1.33(4) 

3.04(15) 

0.061(3) 

4.1/12.0 

895 


Next, we numerically compute the sign function by 
using the minmax approximation by Chebyshev polyno¬ 
mials of Ref. . We solve eigenvalue problems 

by using the Arnoldi method (subroutines of ARPACK), 
and save 0(80) pairs of the low-lying eigenvalues and 
eigenvectors of the Overlap Dirac operator. 

If n_|_ is the number of exact zero modes of plus chirality 
and ri- of minus chirality in the spectrum the massless 
Overlap Dirac operator, the topological charge is defined 
as Q = n+ — n_. The average square of the topological 
charge (Q^), and the topological susceptibility {Q'^) jV 
are computed from the topological charges. All of our 
results are listed in Table im 


C. Eigenvalues and Spectral density 

The eigenvalues A of the Overlap Dirac operator lie in 
the complex plane on a circle of centre (p, 0) and radius 
p as shown in Figure [TJ In our case, p = 1.4. We con¬ 
sider instead the eigenvalues Ximp of the improved mass¬ 
less Overlap Dirac operator D’'^P{p) [2^. They lie on 


the imaginary axis as shown in Figure [T] The improved 
massless Overlap Dirac operator D^™P(p) is defined as: 

'd(p) (II) 

The spectral density p(A, V) is defined as 

P(A, R) = 1(^<5(A-A)). (12) 

A = Im(Aimp)- We show a typical spectral density 
p(A, V) of the non zero modes in Figure [H 


D. The number of zero modes, the topological 
charge, and the topological susceptibility 

In our simulations, we never observe zero modes of 
opposite chirality in the same configuration. The zero 
modes have either all -I- chirality or all - chirality in 
each configuration: as a consequence the number A^obs 
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V = 18^ (3 = 6.00 




FIG. 1. An example of distribution of eigenvalues A, and 
improved eigenvalues Ximp- In this configuration there is one 
zero mode. The right figure is an enlargement of the region 
Ass 0. 



FIG. 4. The average square of the topological charges (Q^) 
versus the physical volume VjrQ. The linear line is (Q^) = 
A ■ V/rl + B, A = 6.8(2) x 10-^ B = -0.20(13). 




FIG. 3. The number of observed zero modes Nzsro versus 
the physical volume V/r%. The continuous curve is Nzero = 
y/A • + B, A = 4.9(3), B = -0.18(5). 


of observed zero modes is Nobs = \n+ — n_| = IQI and 
^obs = Q^- By definition Q = — ri-. 

The average Nzero = {Nobs) is plotted as a func¬ 
tion of the volume V in Figure [3] and is proportional 
to \/V : fitting a function Nzero = \/A • V/r^ -|- i? to 
the data gives A = 4.9(3) x 10“^, B = —0.18(5), and 
X^/ndf = 15.1/15.0. The fitting range in physical volume 
units V/tq is from 24 to 330. (Q^) is consequently linear 
in V [Figure 0 ]. The slope is the topological susceptibil¬ 
ity which is determined to be {Q‘^)rQ/V = 6.8(2) x 10“^. 
This value is consistent with the determination from the 
widths of the distributions of topological charges and 
with the determinations of other groups, as discussed be¬ 
low. 

However our counting of zero modes contrasts with 
the expectation that the number of zero modes be pro¬ 
portional to the volume. Indeed, if zero modes are asso¬ 
ciated to localized configurations like the instantons, and 
the vacuum is made of independent regions of size of the 
order of the correlation length, the number of instantons 
is proportional to the volume. In addition (n+) = (n_) 
by invariance under CP. 

If these arguments are correct the only way out is that 
pairs of instantons of opposite chirality somehow escape 
detection by our way of counting zero modes, at least at 
the rather modest volumes of our simulations. This ef¬ 
fect preserves anyhow the value of the topological charge 
and would not affect the determination of the topological 
susceptibility. We shall use this assumption in the next 
Section to extract from our data the number density of 
instantons. A careful study of this phenomenon, however, 
should be done, to put on safer basis the determination of 
the topological susceptibility. A complementary strategy 
could be to improve numerically the determination based 
on the bosonic sector of the theory [2^, which does not 
rely on the counting of zero modes, and cross-check the 
two methods. 
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FIG. 5. A distribution of topological charges Q. The lattice 
is F = 16"*, ,5 = 6.00. 


FIG. 6. The topological susceptibility at (3 = 6.00 

compared to results of other groups 1^ . 


As anticipated above the topological susceptibility can 
be extracted from the distribution of topological charges. 
A typical distribution of the top olog ical charges is that 
in Figure [5J a Gaussian function 


Pq = 


e 

72^py 


(13) 


fits the data with x^/^df = 7.0/14.0, and topological 
susceptibility {Q'^)rQ/V = 6.7(5) x 10“^. This value is 
consistent with the value = 6.8(5) x 10“^ di¬ 

rectly computed from the number of zero modes. More¬ 
over, seventeen distributions of the topological charges 
are computed from our seventeen different lattices, and 
we fit the Gauss function to all distributions. All the re¬ 
sulting values of x^/nd/ are in the range from 0.3 to 1.3 
[Table [IT] • The topological charges all have the Gaussian 
distribution. 

In Figure m we compare the topological susceptibility 
at the same /3 = 6.00 with results of other groups. The 
figure shows that our results are consistent with them. 
Moreover, we check that the finite lattice volume does 
not affect the topological susceptibility up to L = 2.1 [fm] 
as indicated in Figure [71 


E. Topological susceptibility in the continuum limit 

Last, we fix the physical volume at the value F/rp = 
50.00 marked in Figure [71 and extrapolate the five data 
points of the topological susceptibility to the continuum 
limit using a linear expression (Q^)?'q/F = cq + cio^. 
The results are consistent with other groups. We are 
thus confident that eigenvalues and eigenvectors of over- 



FIG. 7. The topological susceptibilities (Q^)ro/F in our sim¬ 
ulations as shown in Table [III versus the physical volume VfrQ. 
A dotted line marks the physical volume Vfr* = 50.00. 


lap fermions in our simulations are properly computed. 


Our result: x = (1-86(6) x 10^ [MeV])^ 
X = (1-91(5) X 10^ [MeV])'^ 


Ref 

Ref 


28|: X = (1-88(12)(5) x IR^^eV])^ 


(14) 

(15) 

(16) 


The theoretical expectation 29|, [30j : 


X = “ 2m^)|ea;p ^ (1-80 X 10^ [MeV])" 


(17) 


III. INSTANTONS 

We want to compute the number of instantons from 
the number of zero modes, but, as anticipated, there are 
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problems to determine it since some of them escape detec¬ 
tion. In any translation invariant model indeed, e.g. the 
instanton liquid model, the number of instantons linearly 
increases with the physical volume. The number of our 
zero modes instead clearly increases as the square root 
of the physical volume [Figure [3]. We never observe zero 
modes n+ and ri- of opposite chirality in the same config¬ 
uration: our number of zero modes always coincides with 
the topological charge. The distributions of the topolog¬ 
ical charges determined in this way are Gaussian with 
0.3 < x^/ndf <1.3, and agree with other groups results. 

All that suggests that we observe the net number of 
zero modes. At least at our rather small physical lattice 
volumes for some reason pairs of zero modes of opposite 
sign seem to appear as non-zero modes. This can explain 
why we obtain the correct topological charge anyway. To 
estimate the density of instantons we can use a model 
based on the reasonable assumption that the instantons 
of both Chiralities are uniformly distributed in space- 
time and independent. 


A. The number of instantons 


When Q 1, Ni 1 a.t ^ fixed 


P(Q) 


_ q" 
e ^ 


(23) 


Finally, the number of instantons is determined as Nj = 

(Q^) = {NLro)- 


B. The density of instantons 

We obtain the instanton density by fitting a linear 
function Nj = 2pirQ ■ V/tq + B to {Q^) as shown in 
Figure [H All of data points are included in the fit and 
X^/ndf = 11.9/15.0. The slope is 2pirQ = 6.8(2) x 10“^ 
and the intercept is i? = —0.20(13). The intercept is 
compatible with zero. Finally, the instanton density is 
evaluated as 

Pi = 8.3(3) X lO-'^ [GeV^]. (24) 

This result is consistent with the instanton liquid model. 
Ref. 0. 


Let us denote the number of instantons of positive chi¬ 
rality in a volume V by , the number of instantons of 
negative chirality by n_. Of course 

(n+) = (n_) = ^ = PiV, (18) 


Pi is the density of instantons. Because of CP invariance 
instantons and anti-instantons have the same distribu¬ 
tion. The instantons are independent and the distribu¬ 
tion Poisson-like O (): indeed we can view a space-time 
of volume Vst as covered by non overlapping hypercubes 
of size V, which are all equivalent by translation invari¬ 
ance, independent from each other if their size is bigger 
than the correlation length, and can all be used at the 
same time to determine the distribution, which is then 
Poisson-like. 


P(^+)= ^ 

77, I 



(19) 


p(«-) = ^ 


1 / A/ \ " zIIl 


e 2 


( 20 ) 


The resulting distribution for Q = n+ — n_ is 

OO 

P(g) = ^ P(n_)P(n_ +Q) = eM-Ni)lQ{Ni). 

n _—0 

( 21 ) 

Here Iq{x) is the modified Bessel function 


IV. THE MONOPOLE CREATION OPERATOR 

In this section, we briefly recall the definition of 
the monopole creation operator and the method to 
count monopoles. In order to understand the relation 
between instantons and monopoles we add monopole- 
antimonopole pairs of opposite charges in the configu¬ 
rations using the monopole creation operator, and mea¬ 
sure the variation of the number of instantons. We first 
check that the monopoles are successfully added in the 
configurations by counting the additional monopoles. 

The monopole creation operator is defined in 0-0. 
Specifically, in this study, we use the monopole cre¬ 
ation operator defined in 0, [Eq. (41) et seq.], which, 
acting on the vacuum state, produces a static pair of 
monopole-anti-monopole of opposite charge propagating 
from t = —OO to t = -|-oo. The operator is defined as 

p. = exp(—/3AS'). (25) 

AS” is defined by modifying the normal action at the time 
t by replacing the usual plaquette by Il^^(n), as 

follows: 

S + ^= Y, Re(l-n^,(n)) (26) 

= Il^i,(n) in all sites n with no ^ t; at no = t the 
space-space components i,j = 1 — 3 are again unmod¬ 
ified Ilij{t,n) = while Ilio{t,n) is a modified 

plaquette with inserted matrices Mi{n) and Mi{ny, 


k=0 ^ ^ ^ 


( 22 ) 


n,o(t,n) = —^TT[U,(t,n)M}{n + i) 

Tr[/] 

X Uo{t, n + i)Mi{n + i)Ul{t + 1, n)UQ{t, n)] (27) 
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TABLE III. The definition of Distance N and the location of 


the monopole and the anti-monopole in the four-dimensional 
space {x,y,z,t). V = xT. 


N 

Monopole 

Anti-monopole 

Odd 

/L+iV+2 L+N+2 L+1 T\ 

\ 2 ’ 2 ’ 2 ’ 2 / 

/L-N+2 L-N+2 L-1 T\ 

\ 2 ’ 2 ’ 2 ’ 2 / 

Even 

/L+iV+l L+N+1 L+1 T\ 

\ 2 ’ 2 ’ 2 ’ 2 / 

/L-N+1 L-N+1 L-1 T\ 

\ 2 ’ 2 ’ 2 ’ 2 / 


Tr[/] is the trace of the identity, and the matrix Mi{n) is 
the discretised version of the classical field configuration 
A°(n — Xk) produced at site n by the monopoles to 
be added in the locations Xk, namely 

Mi(n) = exp(i^ A°(n - xl)), (i = x, y, z). (28) 

k 


The form used for the monopole fields in a spherical co¬ 
ordinate system (r, 0, /p) centred at the monopole is 
Wu-Yang: 


(z) - z ^ 0 



/ me sin 0(l-|-cos 8 ) 


2gr 


Sin 9 
cos (/>(l+cos 9) 


2gr 


sin 9 

0 


^3 

^3 


[a) Uz — z < 



/ _ rur. sin (^(1—cos 9) \ 

2gr sin 9 ^ 

rric cos (^(1—cos -v 
2gT 's\n 3 




sin 9 
0 


The electric charge is 


9 = 



(gauge coupling constant) 


(29) 


(30) 


(31) 


We give the monopoles magnetic charges 
rric — 0,1, 2, 3, 4. One monopole has charge -|-mc and 
the other has charge —rric- The total magnetic charge 
is zero, rric = 0 is the reference configuration with no 
monopoles added. 

The monopole of charge +mc and anti-monopole of 
charge —rric are placed at time slice t = T/2: the choice 
is irrelevant since boundary conditions are periodic in 
time. 

The locations of the monopoles in the lattice V = x 
T are chosen in the (x, y, z, t) space as in Table Hill where 
also the distance N between the monopoles is defined. 

While Monte Carlo simulations are carried out, the 
pair of monopoles makes long monopole loops in the con¬ 
figurations. 


V. DETECTING THE ADDITIONAL 
MONOPOLES 

To verify that the monopoles are successfully added to 
the configurations, we detect the monopoles in the con¬ 
figurations, as done in Ref. [iS[ni,i3ii. We generate the 


configurations varying the number of added monopole 
charges from zero to four. The distance between the 
monopole and the anti-monopole is defined in Table HUl 
Next, the configurations are iteratively transformed to 
the Maximally Abelian gauge using the Simulated An¬ 
nealing algorithm. To remove the effects of the Gribov 
copies, 20 iterations are carried out in our simulations. 

Abelian link variables, holding 17(1) x 17(1) symme¬ 
try, are derived by Abelian projection from non-Abelian 
link variables. The monopole current is defined for each 
colour direction as follows: 

k’'^,{Jl) = + v) (32) 

The colour index can be i = 1, 2, and 3. rip^(n + v) is 
the Dirac string [3^ . and the monopole current satisfies 
the conservation law, 

= + (33) 

i i 

The monopole density is defined as follows: 

Pmr^ = ^Y^^\klin)\rl (34) 

i n,/i. 

First, we measure the monopole density as shown in Fig- 
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FIG. 8. The monopole density pmro versus the distance 
between the monopole and anti-monopole. The lattice is 
V 13 = 6.00. 


ure [51 varying the distance between the monopole and 
anti-monopole from one to eight, moreover, varying the 
monopole charges rric from zero to four. The monopole 
densities go to plateaus when the monopole and anti- 
monopole are adequately separated. On this plateau we 
fix the distance D between the monopole and the anti- 
monopole, D = 6 and D=8, and measure the monopole 
density. The monopole density does not depend on the 
distance between them as indicated in Table HVl The 























TABLE IV. The monopole density Pm»"o. D6 and D8 indicate 
Distance = 6 and 8 respectively. Nconf. indicate the number 
of configurations each distance. 


p 

V 

rric 

PmCo (D6) 

PmCo (D8) 

TV -P 

6.00 

14* 

0 

0.893(12) 

0.893(12) 

80 



1 

1.026(12) 

1.032(13) 

80 



2 

1.658(14) 

1.660(13) 

80 



3 

2.281(15) 

2.247(17) 

80 



4 

2.991(14) 

3.024(15) 

80 

6.00 

18* 

0 

0.917(7) 

0.916(9) 

80 



1 

0.981(7) 

0.993(8) 

80 



2 

1.307(7) 

1.305(8) 

80 



3 

1.619(8) 

1.621(9) 

80 



4 

1.998(9) 

2.056(8) 

80 


3.5 
3 

2.5 

2 
E 

^ 1.5 
1 

0.5 
0 

FIG. 9. The monopole density pm^o versus the magnetic 
charge rric. Distance = 8. The data points of the normal 
configurations are slightly shifted horizontally from rric = 0 
to distinguish them. 


□ Normal Conf., V = M'' 

■ 

A Normal Conf., V = 18'' 
A V = 18 '' 




monopole densities computed from normal configurations 
of 13 = 6.00, V = 14^ and 18^ are = 0.940(11) 

and PmrQ = 0.919(7) respectively. The monopole den¬ 
sity increases with monopole charges as shown in Figure 
E However, if the physical volume becomes larger, the 
slope becomes lower. As the physical volume becomes 
larger, the number density of the monopoles is rarefied. 
We add the number of monopoles in the configurations 
by increasing the monopole charges. 

Next, we measure the lattice spacing in the same way 
as in Section fll Al The results are listed in Table El and 
show that the additional monopoles do not change it ap¬ 
preciably. 


A. The monopole clusters 

The monopoles are known to form two clusters [TlL 1^ . 
[S^ in MA gauge. The small (ultraviolet) clusters are 



FIG. 10. The histogram of the length of the monopole loops 
L/a. The normal configurations with two flavors of dynamical 
Wilson fermions, V/ro = 1.565(13) x 10®, ro/a = 6.05(5), are 
used. 



FIG. 11. The histogram of the length of the monopole loops 
L/a. Monopoles with charges me ranging from 0 to 4 are 
added to the normal configurations. The lattice is V = 14'*, 
P = 6.00 


made of the short monopole loops, which are short-range 
fluctuations. The large (infrared) clusters, which perco¬ 
late through the lattice and wrap around the boundaries, 
are made of the longest monopole loop in each color direc¬ 
tion. The way how to compute numerically the monopole 
world line in four dimension is explained in [3^ . If the 
physical lattice volume is large enough, the small clus¬ 
ters and the large clusters are separated as in Figure 
nni In quenched SU(2) study the long monopole 
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TABLE V. Determination of the lattice spacing. The lattice spacing and are computed by our simulations, n is the 
number of smearing steps and a is the weight factor of smearing. 


P V 

Distance 

rric 

[fm] 

[fm] 

(n, a) 

FR {Ri/a) 

X^/ndf 

^conf. 

6.00 IS'^ 

6 

0 

0.092(3) 

0.0967(13) 

(30, 0.5) 

1.8 - 8.0 

4.0/4.0 

500 



1 

0.0951(8) 

0.1005(7) 

(25, 0.5) 

0.9 - 8.0 

5.2/5.0 

500 



2 

0.097(3) 

0.1021(17) 

(20, 0.5) 

1.8 - 8.0 

3.8/4.0 

500 



3 

0.097(4) 

0.101(2) 

(25, 0.5) 

1.8 - 7.0 

2.9/3.0 

500 



4 

0.1010(11) 

0.1068(10) 

(25, 0.5) 

0.9 - 6.0 

3.2/3.0 

500 


8 

0 

0.092(3) 

0.0968(14) 

(30, 0.5) 

1.8 - 8.0 

4.0/4.0 

500 



1 

0.095(4) 

0.100(2) 

(25, 0.5) 

1.8 - 6.0 

2.0/2.0 

500 



2 

0.0980(8) 

0.1036(7) 

(25, 0.5) 

0.9 - 9.0 

6.1/6.0 

500 



3 

0.099(4) 

0.104(2) 

(15, 0.5) 

1.8 - 7.0 

3.0/3.0 

500 



4 

0.1009(10) 

0.1062(9) 

(20, 0.5) 

0.9 - 8.0 

3.3/5.0 

500 


loops only exist in confinement phase, and dominate the 
string tension. The long monopole loops are therefore 
considered to play the important role of producing color 
confinement. Going to the Maximally Abelian gauge is 
essential to divide monopole into clusters. 

We create a histogram of the length of the monopole 
loops when one pair of monopoles with charges from zero 
to four are added as shown in Figure [TT] To clarify which 
cluster increases with the monopole charges, we deduct 
the sum of the longest loops from the sum of all loops. 
We define the remainder of the subtraction as the sum 
of short loops. The averages of the sums of the long 
loops and short loops divided by the total number of 
configurations are determined. The results are shown in 
Figure [TH We conclude that the monopoles added are 
long monopoles wrapping the lattice, as expected. The 
length of the long loops is proportional to the charge: a 
charge rric is equivalent to rric monopole pair of charge 1. 


8000 


6000 

“'4000 


2000 


0 


FIG. 12. The average of the long and short monopole loops 
divided by the number of configurations. The lattice is F = 
14"^, P = 6.00. The data points of the normal conhgurations 
are slightly shifted from me = 0 for graphical reasons. 
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VI. THE RELATIONS BETWEEN ZERO 
MODES, INSTANTONS, AND MONOPOLES 

A. Simulation details 


2.5 
2 

§ 1.5 

N 

z 

1 

0.5 
0 

FIG. 13. The number of zero modes Nzsro versus monopole 
charges rric- The data point computed from the normal con¬ 
figurations is slightly shifted from rric = 0. 


■ Normal Conf. 
□ Distance = 8 
A Prediction 




m. 


We generate configurations with one monopole-anti- 
monopole pair added with different magnetic charges rric- 
The distances between the monopole and anti-monopole 
are fixed at 6, and 8: slightly changing the distance be¬ 
tween them, allows to check finite lattice volume effects. 
The simulation parameters and the data are listed in Ta¬ 
ble ED The Overlap Dirac operator is constructed from 
the gauge links of the configurations. The eigenvalue 
problems are solved by the Arnold! method, and 0(60) 
pairs of the low-lying eigenvalues and eigenvectors of the 
Overlap Dirac operator are saved. We then count the 
number of zero modes, and calculate the average square 
of topological charges, i.e. the number of instantons. Fi¬ 
nally, we compare the simulation results with the predic¬ 
tions based on the hypothesis that the added monopoles 
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FIG. 14. The average square of topological charges (Q^) ver¬ 
sus monopole charges rric. The data point computed from the 
normal configurations is slightly shifted from rUc = 0. 
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TABLE VI. Charges of the added monopole, together with 
number of zero modes Nzero, and the average square of topo¬ 
logical charges (Q^). The lattice is /3 = 6.00, V = 14“*. 


Distance 

rric 

NZero 

(Q^) 

^conf. 

6 

1 

1.50(6) 

3.6(3) 

400 


2 

1.78(7) 

5.4(4) 

401 


3 

2.03(8) 

6.7(5) 

423 


4 

2 .10(8) 

7.2(5) 

440 

8 

0 

1.39(5) 

3.2(2) 

501 


1 

1.55(6) 

4.0(3) 

502 


2 

1.85(7) 

5.6(4) 

510 


3 

2.09(7) 

6.8(4) 

500 


4 

2.16(8) 

7.5(5) 

505 


do not perturb the vacuum and specifically the distribu¬ 
tion of instantons, but can only change the total topo¬ 
logical charge. 

Note: We do not do smearing, cooling, or MA gauge 
fixing in simulations. The number of zero modes, Dis¬ 
tance 6 rric = 0, coincide within errors with Distance 8. 
We take this as a check of volume independence. 

In Figure [TS] and Figure [T3] we show the number of 
zero modes and the average square of the topological 
charges respectively as functions of the monopole charge 
rric- rric = 0 is the ordinary case with no monopoles 
added. The results are compared with the predictions 
developed below. 


B. Predictions 

The creation operator Eq. (1251) acting on the vacuum 
produces a pair of static monopoles propagating in time 
from — oo to -Poo. If the vacuum is a dual superconductor 
it is a coherent state and an eigenstate of the creation op¬ 


erator of the pair, and is therefore left unchanged. In any 
case, as a dual superconductor the vacuum shields the 
monopoles, so that they are one-dimensional structures 
and do not influence the external space time: 0{V~'i). 
For these reasons we expect that the distribution of zero 
modes Eq. (HU), (HOI) is unaffected by the additional 
monopoles especially at large space-time volumes V, ex¬ 
cept for the addition of some instantons. 

If each monopole-antimonopole pair of charge rric = 1 
produces q instantons or anti-instantons with equal prob¬ 
ability (CP invariance), and the original distribution is 
not modified, the value of (Q^) becomes 

((Q + AQ)") = (Q2) + 2{Q){AQ)i + {(AQf)i (35) 

Here (• • •) means average on the configuration, (■■■)/ 
the average on the charge distribution of the additional 
M = rric X q instantons. The middle term in the right- 
hand side of Eq. (1551) vanishes since (Q) = 0. To define 
(•••)/ we assume that each chiral charge can be ±1 with 
equal probability. The probability to have k positive and 
M — k negative instantons is 

P(k,M) = ^{^^y (k = Q...M). (36) 

Corresponding to this partition AQ — 2k — M so that 

M 

{(AQf)i = '^P(k,M)(2k-Mf = M. (37) 

0 

The sum is easily computed as shown in Appendix 
Recalling that M = rricq, the prediction follows that, in 
presence of rric pairs of monopoles 

{{Q + AQ)^) = (Q^) -\- rricq (38) 

Comparing with the data [Figure[T4] it follows that 9 = 1. 
Fit of a linear function (Q^) = A-rric + B to the results of 
the average square of the topological charges computed 
from simulations gives a slope A compatible with 1, and 
intercept B compatible with (Q^) = 3.17(19) (the value 
from normal configurations) [Table IVIII| . 

The observed behaviour supports the assumption that 
monopoles do not alter the vacuum. Each pair of 
monopoles adds one instanton or one anti-instanton with 
equal probability. Having fixed once and for all that 9=1 
we can explicitly compute the topological charge distri¬ 
bution, under our basic assumption, and compare with 
data. Since our monopole configuration is CP invariant, 
the distribution has to be even in Q. In the simple case 
rric = 1 the distribution is the sum with equal weight ^ of 
two ’’unperturbed” distributions Po(Q) centred at Q = I 
and Q = -1, P{Q) = \[P^{Q + 1) + Po(Q - !)]• Pn(Q) 
is well approximated by Eq. (1551) . In the case of general 
rric 

P(Q,mc) = 

1 / \ 

= T.[k) + + PoiQ -mc + 2k)] 

(39) 
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TABLE VII. The results by fitting the function Eq. (1391) . 



rric 

Ni 

X^/ndf 

^conf. 

Distance 6 

1 

2.5(3) 

8.5/11.0 

400 


2 

3.3(4) 

9.9/14.0 

401 


3 

3.3(4) 

15.5/15.0 

423 


4 

3.0(5) 

8.5/16.0 

440 

Distance 8 

1 

2.7(2) 

13.6/13.0 

502 


2 

3.4(4) 

8.2/15.0 

510 


3 

3.6(4) 

14.1/14.0 

500 


4 

3.1(5) 

16.2/17.0 

505 

TABLE VIII. The final results. 


A 

B 

FR (me) 

X^/ndf 

Prediction 

1 .00(6) 

3.2(2) 

1 - 5 

0.0/3.0 

Distance 6 

1.05(11) 

3.06(18) 

0 - 4 

5.7/3.0 

Distance 8 

1.13(10) 

3.15(18) 

0 - 4 

2.2/3.0 


The distribution is CP invariant if Pq is such. 

By use of Eq. (IMl) we can predict (|Q|) numerically, 
using the topological charges computed from the normal 
configurations as input value, keeping in mind that pairs 
of zero modes with opposite chirality escape detection. 
The computation is done in detail in Appendix section IbI 
and compared to the data in Figure [131 

The histograms of the topological charges Q for each 
monopole charge rric are compared to the prediction ob¬ 
tained by multiplying the distribution Eq. (IMll x the 
number of configurations in Figure 1151 The input value 
Ni = 3.17 is used. The topological charge distribu¬ 
tions computed by simulations are represented by the 
histogram, and the prediction by the continuous line with 
symbols. The y'^/ndf indicates agreement with the pre¬ 
diction by Eq. dMl)- 
In the general case 

{Q^) = {S^)+m, (40) 

((5^) is nothing but (Q^) at rUc = 0. Details of the com¬ 
putations are given in Appendix section iBl 


VII. SUMMARY AND CONCLUSIONS 

The aim of this paper is to shed light on the relation 
between instantons and monopoles. 

Our system is an S'[7(3) pure gauge theory on the lat¬ 
tice. Our strategy is to add pairs of static monopoles of 
opposite magnetic charge to the system and detect the 
consequent variation in the number of instantons. 

To add monopoles we use the creation operator devel¬ 
oped in Ref. If tlie vacuum is a dual superconduc¬ 

tor we expect that adding the monopoles does not modify 
the ground state, except in the number of monopoles. We 
check by the technique of Ref. Bini that we are really 



FIG. 15. The comparisons of the distributions of the topolog¬ 
ical charges computed by simulations (Distance = 8) and the 
prediction. The monopole charges rric are from 1 to 4. 


adding monopoles: we find that we are indeed adding 
’’long” monopoles, wrapping around the lattice in the 
(periodic) time direction. 

To detect instantons we count the zero modes of a 
fermion overlap operator. In doing that we realise that 
this tool is only able to detect the net topological charge: 
the zero modes in each configuration are either all of posi¬ 
tive chirality, or all of negative chirality, but never appear 
with opposite sign. The number of observed zero modes 
is proportional to the square root of the space-time vol¬ 
ume and (Q^) to Vst as expected. The only way to 

reconcile this fact with the obvious expectation that the 
number of zero modes be proportional to Vst, the vacuum 
being translation invariant and the correlation length fi¬ 
nite, is to assume that, for some reason, pairs of instan¬ 
tons of opposite chirality escape detection. Maybe this 
phenomenon disappears in larger lattices. This should be 
anyhow quantitatively understood, also to have on safe 
ground the determination of the topological susceptibil¬ 
ity. Based on this model we are able to determine the 
number of instantons and its variation due to the addi¬ 
tion of the monopoles. 

We find that an instanton or an anti-instanton are 
added to the system with equal probability by each pair 
of monopole-anti-monopole of magnetic charge 1. The 
variation of the chiral charge is equal to the magnetic 
charge of the monopole of the pair. This is the our main 
result. 

We keep running simulations in order to further clar¬ 
ify that Chiral symmetry breaking is induced by the 
monopole-anti-monopole pair Ref. |13j . 
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Appendix A: Proof of Eq. (1381) 


Consider the function f(x) defined as 

/(x) =^P(A,M)x'=(i)“-'= = (Al) 

k 


Compute 5 ( 1 ) with g{x) = on both sides of 

the last equality. The result is Eq. (1381) . 


Appendix B: The computation of (|Q|) 


Erom the distribution of topological charge Eq. dMl) we 
can compute (|QI) as a function of rric- |Q| is the observed 
number of zero modes, being our detection blind to pairs 
of instantons of opposite chirality. We already know from 
Eq. (1381) that {Q'^) = (d^) +mc where the average (d^) is 
weighted by Po(<5)) the ’’unperturbed” distribution. (<5^) 
is the value at rric = 0 . 

Eor rric = 1 we have 


= ^{(l^ + l|) + (l^-l|)}- (Bl) 


Now, 

(l^ + l|) = 

1 




r*-l 


{5 + l)e~^d5 - {S-\- l)e~^d6 

) J-i 


g 27^ — _ 

TT V2ttN 


r-1 


e dS — 


~^dS 


'-1 


U-oc 

fm 1 


(B2) 


{5 - l)e~^d6 - J {6-l)e~^d5 

r*l o poo o 


1 




e 2 NdS - 


e 

(B3) 


The number of zero modes Nzero that we observe in our 
simulation is 

1 


Nzero =^{{\S+l\) + {\5-l\)} 

'm 


V^ttN J-i 


'^dS 


(B4) 


The integrals are computed numerically, by substituting 
N = (d^) = 3.17(19) of Nconf. = 592 computed by our 
lattice simulations. The result is Nzero = 1.64(6). 

For rric = 2; 


Nzero — ^{(115 + 21 ) -|-(|(5-2|)}-|-- 


1 


V2nN 
= 1.82(8) 


Ne ^ + J e ^^d6} + 


N 


V2TrN 

(B5) 


For rric = 3; 

iVze.o = i{(|5 + 3|) + (|d-3|)} + ^{(|<5 + l|) + (|d-l|)} 


\/2'nN J -3 


1 r 4V _JL 6 

= —< - P 2N -I- 

8 1^2^ 

3 r 47V 2 

H—< -e -|- 

8 i V2^ V2^ . 

= 1.99(9) 


e ^^dS 


'^^d6 


(B 6 ) 


For rric = 4; 

Nzero = A 

+ ^{(l^ + 2 |) + (l^- 2 |)} + ^(|d|) 


1 , 

f 4V , 8 

j 

“ 16 ' 

1 V27rV 

V2ttN . 

+ 4] 

' 4V 2 4 

I 

^ e~^d5^ 

^y/2-KN 

y/2TTN J 


= 2.15(9) 


3N 


NZero 


(B7) 
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